November 22, 2012 4;14 Quantitative Finance paper 

Quantitative Finance, Vol. 00, No. 00, Month 200x, 1-18 



Impact of meta-order in the Minority Game 



o 

(N 

> 
O 

(N 



A C Baratof, I Mastromatteo|, M Bardoscia*t and M Marsilif 

t The Abdus Salam International Centre for Theoretical Physics, Trieste 34014, Italy 
I International School for Advanced Studies, Trieste 34136, Italy 

(Received 00 Month 200x; in final form 00 Month 200x) 



H 
d 

(N 

> 
oo 

O 



Wc study the market impact of a meta-order in the framework of the Minority Game. This 
amounts to studying the response of the market when introducing a trader who buys or seUs 
a fixed amount h for a finite time T. This perturbation introduces statistical arbitrages that 
traders exploit by adapting their trading strategies. The market impact depends on the na- 
ture of the stationary state: we find that the permanent impact is zero in the unpredictable 
(information efficient) phase, while in the predictable phase it is non-zero and grows linearly 
with the size of the meta-order. This establishes a quantitative link between information effi- 
ciency and trading efficiency (i.e. market impact). By using statistical mechanics methods for 
disordered systems, we are able to fully characterize the response in the predictable phase, to 
relate execution cost to response functions and obtain exact results for the permanent impact. 

Keywords: Minority Game; Market Impact; Meta-order; Market Microstructure. 



Introduction 



X: 



Market impact is the effect on the price caused by an order, and it measures th e tendency 
of pr ices to move up (down) subsequent to a buy (sell) trade of a given size (jBouchaud 
20l3). Understanding market impact is a fundamental problem in fi nancial markets, and 
it ha s recently been the subject of severa l stud ie s, both ernpirical (IHasbrouck and Seppi 



2010l. IPlerou et al\ l2002l. Ihvons and EvansI l2002l . Ihillo et a.l 



Farmer et al\ l2005l . lOabaix al\ l2005l . iBouchaud et a.l\ l2008l . 



Moro et aLll2009l. iToth. Eisler adl201ll) and theoretica l (lTorre|ll997l. ICrinold and Kahn 



200l lA lmgre n et al 



Toth. Lemperiere et al. 



2005, 



-2000I . 

201 ih . 



From a theoretical perspective, market impact is relate d to the problem of understand- 
ing how information is incorporated into security prices ( Bouchaud et al. 2008 . Kylel 19851 . 
Glosten and Milfrom|[l985l . lHasbroucklll99ll ). 

For practitioners, controlling market impact is essential in order to limit execution costs, which 
arise because each trade impacts adversely on the price. In addition, traders often face the 
problem of executing orders whose volume is much larger than the outstanding liquidity which 
is available for trading at each moment in time, which is typically a factor 10~^ smaller than the 
volumes exchanged daily ( Bouchaud et al. 20081 ). A common strategy consists in splitting large 
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orders in a sequence of child orders that are executed incrementally. This sequence of trades is 
commonly called meta- order. The trading activity generated by meta-orders has been argued 
to be one of the origin s of long range auto-correlation in the flow of orders in financial markets 
( Bouchaud et alJ 20081 ). and generates market impact which extends beyond that of single orders. 
One robust empirical finding is the fact that ma rket impact is a concave fun ction of the order 
size both for individual orde rs ( Bouchaud et al\ 2008 ) and for me ta-orders |Moro et al. 20091 . 
Toth. Lemperiere et al. 20 111 ), in contrast with early models (e.g. Kvld 19851 ) which predict a 



linear impact. 

Understanding market impact requires addressing the strategic interaction between liquidity 
providers and liquidity takers, and in particular mo deling how t he former exploit the predictable 
trading patte rns of the latter. Ky le's seminal paper (lKvlelll985l l and its most recent extension to 
meta-orders (jFarmer et aLll201ll ) build on this key ingredient, but are limited to simple (repre- 
sentative) agent settings. In addition, the assumption of perfectly informed agents hardly suits 
to realistic setting^. More recent the oretical studie s on market impact have been done within 
the framework of order-book models ( Slanina 20081 ). which provide an accurate description of 
the market microstructure. We take a complementary approach by using a coarse-grained de- 
scription of the market, focusing on the collective behavior of a set of heterogeneous adaptive 
market participants. 

The natural framew' o rk for addressing this issue is the Minority Game (MG) 
( Challet and Zhandll997l . IChallet et ai]l2005l . ICoolenll2005l l. which is a paradigmatic model of 
heterogeneous adaptive agents who repeatedly interact in a market setting using simple "trading 
strategies" to either buy or sell an asset, with the objective of being in the minority side (i.e. sell- 
ing when the majority buys or vice- versa). In spite of its simplicity, the MG offers a picture of a 
financial market as a many-agent interacting system where speculators "predate" on predictable 
patterns introduced by other agents, much as in real financial markets. As such, the MG shows, 
in a simplified setting, how markets aggregate information as a result of speculative activity. In 
brief, one can distinguish an "asymmetric" phase, characterized by the presence of statistical 
arbitrages when speculative pressure is low, from a "symmetric" pha se, when there are many 
speculators and the market is unpredicable (jChallet and Marsilil[l999l ). Close to the transition 
between the two regimes several stylized facts observed in f inancial markets can be successfully 



descr i bed by slightly generalizin g the original MG problem (IChallet et aLll200ll . iBouchaud et al. 
200ll . IChallet and MarsilillioO^ ). 



The MG is the appropriate setting for studying the market impact of meta-orders, because 
it reproduces the response of adaptive traders to perturbations introduced by the persistent 
activity of meta-orders. A meta-order is introduced in the MG by adding an extra agent that 
consistently buys (or sells) for a finite time period, after which the extra agent is removed, and 
the meta-order finishes. In order to analyze how this perturbation affects the price, one needs 
to make an assumpt ion on how th e price is related to excess demand. Here we assume a market 
clearing mechanism ( Marsili 200ll ) taking a linear relation between them. Alternative approaches 
are discussed in Sec. [H This allows us to estimate the impact of the meta-order on price, during 
and after its execution. Our main findings are: i) m. the predictable phase the permanent impact 
is nonzero, while in the unpredictable phase there is no permanent impact; ii) the permanent 
impact can be computed analytically, relating it to response functions. Finally, in) we find that 
the market impact is linear with the size of the order, as a consequence of the assumed linear 
relation between price and excess demand. 

The organization of the paper is as follows. In the next section we define the model, discuss 
its phase transition and explain how the meta-order is introduced. In Sec. [3] we present our 



Typically, liquidity takers are perfectly informed on how much they need to trade, but not on the future price of assets. 
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results and in Sec. |4]we draw the main conclusions. Furthermore, we add appendices with the 
calculations used to obtain the analytical results. 



2. Grand-canonical Minority Game 
2.1. Model definition 

The MG describes a set of agents interacting in a market over a sequence of periods labeled 
by a discrete time variable t S N. Each period corresponds to the trading activity in a time 
window \t,t + At). A realistic value for At must be much larger than inter-trade times, e.g. 
At ^ 10 s, and can be thought as being the typical time over which meaningful information 
accumulates, e.g. few minuteo Speculators act in the market by using trading strategies that 
process the value of a public information signal /x and submit orders accordingly. In practice, 
a number of trading strategies is constantly being evaluated and the most profitable are used. 
Profitable strategies are those that dictate to buy when the majority is selling or vice-versa, i.e. 
that place the trader in the minority side as frequently as possible. This captures in a simple 

manner the market making activity of liquidity providers; 

Formally, we consider the grand-canonical MG setup ( Challet et al. 20051 ). The information 



to which traders react can take P values, labeled by jj. = Every time step a new 

signal /i(t) is drawn at random. We s uppos e there are two types of traders: liquidity takers - 



also called producers in lChallet et al\ ()2005l ) - and speculators or liquidity providers. There are 
Np producers and Ng speculators and their only difference is that while the former react to 
information signals in a deterministic manner, the latter are adaptive. 

Specifically, each agent i is endowed with a trading strategy {af}^=i, that dictates whether 
to buy {a^ = +\) or to sell (a^ = — 1) when public information is /x(t) = /i. These strategies 
are independently drawn for each player, uniformly in the set of all 2^ possible binary functions 
/ : {!,..., -P} —7- {±1}. However, while producers always trade according to their strategy, 
speculators may choose to trade or not, adapting their behavior to the profitability of their 
strategy. This choice is encoded in a dynamical variable 4>i{t), which takes the value 4>i{t) = 1 if 
agent i decides to trade, and (j)i{t) = otherwise. The profitability of a strategy is quantified by 
a score Ui{t) that is updated according to 

Ui{t + l) = U,{t)-af^A{t), (1) 
where the excess demand A{t) is given by 



A(t) = ^f) + J]afV.(t) (2) 

i=l 



and Ap is the sum of the contributions from producers when /x(t) = /x. The decision of playing 
or not is taken in accordance to the score of the speculator. In the simplest case, agent i does not 
trade if Ui{t) < whereas she trades if the score is positiv^. In summary, a time-step iteration 



^One method to estimate time-scales over which i nformation flows is that of studying the dynamics of cross-correlations 
between different stock price s. On very liquid stock s, iMastromatt co et oil ||2011| ) show that for At < 1 min price dynamics is 
dominated by noise, whereas iBorghesi et al\ l|2007f ) show that the structure of correlations of dai ly returns can be r ecovered 
already from the analysis of returns over At = 5 minutes. Similar results have been d iscussed in lToth et al.\ | |2006| ). 
■^This rule can be generalized using a stochastic choice model llChallet et aZ.I [20051 ). without changing the main results 
presented in this paper. 
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of the model is realized in the following way: 

- randomly generate a information variable G {1, . . . , P}; 

- for each speculator i, set (j)i{t) = 1 if Ui{t) > 0, and (j)i{t) = otherwise; 

- calculate A{t) using Eq. ([2]); 

- for each speculator i, update Ui{t + 1) by Eq. ([TJ; 

- set t — ^ t + 1. 



2.2. Phase transition in MG 

As discussed in IChallet and Marsilil ^200^ ) this model has a very rich phenomenology. In the 
limit of large markets, analog to the thermodynamic limit in physics, the non-trivial regime 
arises when P,Ns,Np — )• oo, while rig = Ng/P and rip = Np/P remain finite. A central role in 
the characterization of the collective behavior is played by the quantity 



(3) 



where 



{A\^^) = A^ + Y,<{<P^ 



(4) 



1=1 



denotes the average of A{t) conditional on fi{t) = fj,. H detects the presence of statistical arbi- 
trages, in that ii H > there is at least one value of fx which allows for a statistical prediction 
of the excess demand when fj,{t) = /u. For this reason H is called predictability. A remark- 
able feature of the MG is that the dynamics converges to states of minimal predictability H 
(jChallet et a?.ll2005l ). This is not only appealing conceptually, as it says that speculators leave 
the market as unpredictable as possible given their strategies, but it also opens the possibility 
for analytic approaches to the MG. Indeed, it turns the study of the stationary state of the 
MG into the analysis of the minima of H, as a function of the variables This is akin to 
studying the ground state properties of a disordered system in physics, a problem for which a 
plethora of powerful tools has been developed in recent years (|Challet et a;.ll2005l . ICoolenll2005l l . 
In Appendix |A] we show how this is done by calculating the critical line and the predictability 
in the asymmetric phase. In brief, one finds that 



H Up + Gug 



(5) 



where G is given in Eq. ()A9p in the Appendix [XI and 



X 



1 ^= 

—y 

Ng ^ 



6 at 



(6) 
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can be interpreted as an agent susceptibility, as it measures the response of agents' behavior 
to a change a-* — >• + 5oj in their strategieql|. The relation between x and G on and Up is 
discussed in the appendix, here we focus on the resulting picture of the behavior of the MG, 
which can be summarized as follows. For a fixed value of n^, if is small enough the market 
is predictable {H > 0) and the susceptibility x is finite. In words, the number of speculators is 
not enough to exploit the amount of statistical arbitrages introduced by producers. When Ug 
increases, H decreases whereas x increases, in agreement with Eq. ([5|). Beyond a critical value 
n* the susceptibility diverges {x = oo) and the market is unpredictable {H = 0). The point n* 
marks a phase transition in that it separates two regions with distinctly different features. In 
this paper we focus on the case Up = 1, for which the critical value is ?i* 4.15. 



2.3. Meta-orders in the Minority Game 

In order to introduce a meta-order in the grand-canonical MG we add a fixed buyer to the system. 
That is, we add one extra producer that always buys independently of the public information 
pattern. This extra producer buys a quantity h of assets at each time, for T time steps. This 
models a meta-order of duration T, and total size hT. Since typical relaxation times in the 
MG are proportional to P, we'll consider T proportional to P in what follows, which might 
correspond to time-scales of one trading day. Note that turning on the meta-order corresponds 
to changing the excess demand to 

Af,{t) = he{T - t)e{t) + + afU^{t) , (7) 

i=l 

where 0{x) = 1 for x > and 0[x) = otherwise. This way of modeling the unwinding of a meta- 
order is reminiscent of a very common scheduling strategy (the TWAP schem e), in which the 



trade r's target is to execute a quantity of shares which grows linearly in time (iBerkowitz et al 
1988h . 



It can be shown that, impo sing a market clearing condition (see iMarsilil 12001 



De Martino and Marsili et al. 20061 ). the quantity A(t) is proportional to the return of the log- 



arithm of the price, which is given by 

logp{t + l) = logp{t) + ^A{t). (8) 

Different market mechanisms can imply a different relation between price change p{t + l)/p{t) 
and excess demand A[t). This relation has implications in how the market impact depends on 
the order size, which is of central importance in market impact studies. Our definition is likely 
to be appropriate at large time scales (say, between one day and a week), in which one expects 
stochasticity not to play a dominant role in the description of the market clearing mechanism. 
We shall return to this point in the conclusions. We define the price impact as the averaged 
log-price change in presence of the meta-order, minus the trend component which would have 



^This is a standard result l IMezard et a^.l [19871 ). see e.g. App endix Al.l in IMarsilil j2009l ) for an explicit derivation in 
an interacting agents' context. In the dynamic theory of MG llCoolenI 120051 ). x emerges as the integrated response of an 
instantaneous perturbation. Loosely speaking, it measrues the time over which agents's behavior maintains memory of 
perturbations. 
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been present in the unperturbed case. This corresponds to the quantity 

^W = ^ E (Ms) - A,=o{s)) , (9) 

0<s<t 



where s = corresponds to the beginning of the meta-order, the average over the stationary state 
distribution is denoted with (• • • ) and the one over the possible choices of the fixed strategie^ 
has been written as (• • • ). This averaging procedure mimics the averages in empirical studies , 
which are on different orders at different times, as wel l as on different markets ( Moro et al. 20091 . 
Toth. Eisler et al. 201ll . lT6th. Lemperiere et al. 2011 ). Notice that [A^^q) is expected to be zero 



only after averaging over the fixed strategies, as in the asymmetric phase its value is finite, and 
fluctuates from realization to realizationo The averages in Eq. ^ are numerically estimated 
according to the following procedure. Initially all the scores are set to zero and the model is run 
until it relaxes to the stationary state. In this initial phase of the evolution speculators adapt 
to the information injected by producers in the absence of a meta-order. Then we compute 
(^/i=o) taking a time average in the stationary stat^. After that the meta-order is turned on 
and we start measuring the impact. Finally the average over the fixed strategies is computed by 
repeating this procedure for a large number of realizations. 

We will also consider A(t) for t > T, since we are interested in the response of the system 
after the meta-order is turned-off. Of particular interest is the permanent impact, which reads 
A* = limt_5.oo and the average execution cost 



(10) 

which measures the mean price increment incurred during the meta-order. In particular, it is a 
simple exercise to show that the cost of an execution for a trader is related to the maximum 
price increment p(0)A(T) as 



p(0)A=p(0) 



t=l 



A(r) {p{T)-pm 

The estimation of the above quantity provides a practical mean to evaluate the concavity of price 
impact: for a market impact of the form A{t) oc t", one finds that A/A(T) = 1/(1 + a), which 
indicates that concave impact functions have to be associated with a value of A/A(T) > 1/2, 
while for convex impact functions a value smaller than 1/2 has to be expected. 



3. Results 

Let us first analyze how the impact changes depending on whether the model is in the symmetric 
or the asymmetric phase. In Fig.[T]we plot the impact A(t), obtained from numerical simulations. 



^Thc strategics arc fixed in the sense that they do not change during the whole evolution of the system, i.e. before the 
beginning of the order, during the order and after its completion. In physical jargon such variables would be called quenched 
strategies. However, fo r the sake of clarity , we will continue to refer to them as fixed strategies. 

^It can be estimated llChallet et aL|[2005D that (A^— q) is finite for any realization of the MG, and it is described by a 
Gaussian random variable of mean zero and variance H/Ns- 

•^The average has to be taken on times much longer than those over which we observe the impact A(f). 
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Figure 1. Comparison between the MG in the asymmetric and symmetric phases for P — 400, T = 5P, 
and 5 x 10'^ realizations. On the left side we have Hs = 1 where the predictability and the permanent 
impact are non-zero. On the right side = 5, where H = 0. 



as a function of time for = 1 (predictable phase) and lis = 5 (unpredictable phase). In 
the predictable phase the impact grows linearly during the meta-order. After the meta-order 
finishes, it decays to a constant faster than a power law, but with a long relaxation time. In the 
unpredictable phase the situation is rather different. During the meta-order the impact grows 
initially very fast and then it reaches a saturation value. After the meta-order is finished the 
impact quickly relaxes back to zero. Therefore, as a main result of this paper, we find that if the 
market is unpredictable {H = 0) the permanent impact is zero and if the market is predictable 
(H > 0) the permanent impact is non-zero. We point out that the permanent impact is zero in 
the symmetric phase as long as the the duration of the meta-order T is big enough such that 
the impact saturates before the end of the meta-order. If T is much smaller than the time for 
saturation to set in, the permanent impact in the symmetric phase is non-zero. 

In terms of the excess demand we have the following picture: immediately after the order is 
turned on, the (average) excess demand jumps from to h. During the execution of the meta- 
order it relaxes to a value smaller than h. In the unpredictable phase, this value is zero, whereas 
in the predictable phase it is finite, implying a linear growth of the impact. At the end of the 
meta-order (t = T) the excess demand suddenly drops by /i, thus becoming negative. Finally it 
relaxes back to zero, regardless of the phase. 

This result can be explained in the following way. If predictability is zero, it means that the 
market is information free and the system is able to absorb the extra information of the extra 
fixed buyer. Therefore, the value of the excess demand relaxes back to zero during the meta- 
order. On the other hand, when the market is predictable the number of speculators is not 
enough to consume all the information available. Hence, when the perturbation is added, the 
system reaches a new stationary state where the value of the excess demand is non-zero. 

Fig. [2] shows that the impact grows linearly with the order size Q = hT. Indeed the rescaled 
impact A/h collapses on the same curve, when plotted against t/T. This is a consequence of the 
assumed market mechan ism and of linear response the ory, as we shall see, and it contrasts with 
results of empirical (e.g. Toth. Lemper iere et al\ 2011 ) which find a square root law A ~ ^/Q. 
We will come back to this issue in the concluding section. 

These results are fully consistent with the theoretical predictions. As we show in Appendix 
[B| the stationary state of the system in the presence of the meta-order can be calculated by 
minimizing the predictability H with the modified excess demand ([7]) . This allows us to calculate 
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Figure 2. Market impact as a function of time for different values of h, for Ug = 1, P = 400, T ~ 5P, 
and 5 x 10'^ realizations. The collapse of the curves shows that the market impact grows linearly with h. 




Figure 3. Comparison between the saturation value of the aggregate in the presence of the mcta-ordcr 
(Ah) obtained analytically (full line) and from numerical simulations (open circles). The error on the 
points from numerical simulations is at most of the size of the open circles. 



the saturation value of the excess demand before the meta-order is finished, which is given by 



(Ah) = {Ah=o) + 



1 + x' 



(12) 



In order to obtain it numerically we calculate the slope of the curves of A(t) when the linear 
growth of the impact sets in, after some initial transient. In Fig. [3] we can see that the numerical 
results are in full agreement with the exact calculation. 

It is worth noticing that {Ah)/h, which measures the response of the market to a perturbation 
A ^ A + h, is inversely proportional to the susceptibility x of agents' behavior. The market 
impact is larger when agents are not very susceptible to perturbations, whereas it is minimal in 
the symmetric phase {H = 0), when agents respond strongly to any perturbation. 

In order to address the dynamic effects of market impact, it is necessary to reformulate the 
dynamics in a continuum time description which acknowledges the fact that characteristic times 
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3- 



2- 




Figure 4. Comparison of numerical simulations with the permanent impact given by formula (jlSp . The 
numerical results were obtained with P = 400, T = 5P, and 5 x 10'^ independent realizations. The straight 
horizontal lines indicate the value of the permanent impact obtained from ()15|) . 



in the MG are of order P (jChallet et al\\20(]^ . ICooleiill2nn,^ ;i. As long as the perturbation h is 



much smaller than A \/P, linear response theory holds, i.e. 

1 



t/P 

{Ahit)) = {Ah=oit)) + I dTK{t/P-T)h{T), (13) 

1 + X 



oo 



where /i(t) = h{PT) = h6{T/P — t)9{t) is the perturbation and K{t) is a kernel that contains 
the relaxation dynamics. Note that if t, T — >• oo we expect to recover the stationary value of Eq. 
(|12p . hence /_ dTK{T) = 1. Also causality requires K[t) = for r < 0. We remark that in 
order to reproduce the small time behavior of the model (i.e. {Ah{t) — Ah=Q(t)) ~ h for small t), 
the function K(t) must have the form 

K{T) = {l + x)S{T)-xKr{r), (14) 

where Kj.{t) is non-singular in zero, causal (X,.(t) = for r < 0) and J°°^ dTKr{T) = 1. Kr{t) 
in Eq. (fH|) captures how adaptive traders detect the statistical arbitrage introduced by the 
meta-order and react to it, by removing it as a result of their activity. If in addition we assume 
that Kr{T) decays fast enough when r — )• oo, then within linear response theory, we show in 
Appendix [O that the permanent impact is 

The assumption on K{s) is fully consistent with simulations, where the excess demand is found 
to relax faster than a power-law, which is expected since we are not considering the model 
at criticality. In Fig. |4] we see agreement between Eq. (jl5p and numerical results. However, 
relaxation times are rather long even for values of considerably far from the critical value 
n* ~ 4.15, making the error bars become very large before A is saturated. 
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Finally, the average execution cost, which is derived in Appendix [Dj is given by 



A(r) 



dxx'^KixT/P). 



(16) 



For a kernel K{t) of the form (jl4|) and execution times T/P much smaller than the characteristic 
relaxation time of Kt.{t), this quantity is well approximated by 



A _ 1 / X 

A{T) ^ 2 + 12P U + X 



KriO), 



(17) 



while at large times we obtain a value A/A(T) = 1/2 which signals a linear impact. This should 



be compared with A/A(r) ~ 0.6 -r- 0.7 reported in empi rical studies (e.g. [Farmer et al.l 12011 
Toth. Lemperiere et al. 201 ll ). and with the arguments in Toth. Lemperiere et al. ( 201ll ). 



4. Conclusion 

The market impact has been the subject of intense discussion recently. Propo sals of universal 
laws have been put forward ( Farmer et al. 2011 . Toth. Lemperiere et al. 201 ll ) on the basis of 



empirical analysis and simple models. Here we address this iss ue within the framework of th e 
Minority Game. At odds with zero- intelligence (or e-intelligence ( Toth. Lemperiere et al. 20 111 )) 



models, the MG retains a level of rationalit y of the agents with out making strong assumptions on 
their information on future returns (as in Farmer et al. 201 ll ). Rather, traders learn adaptively 



on the profitability of their strategies. In addition, the minority rule is rather appropriate if one 
focuses on market making activity, where the profit of agents comes from matching imbalances 
in demand. Finally, the MG is the ideal platform to address how markets digest information on 
order imbalances from meta-orders and to address the issue of the market impact beyond its 
immediate consequences. 

We have shown that the MG provides realistic curves for the market impacts and sheds light 
on the mechanisms which are responsible for it. Most interesting, we derive a relation between 
market impact and information efficiency, whereby the permanent impact of a meta-order is 
smaller the more the market is free of statistical arbitrages. For a perfectly "symmetric" market 
where returns are statistically unpredictable, prices should revert to the values before the order 
started, and hence the permanent impact should be zero. Understanding whether this relation 
extends to other modeling frameworks or whether it is confirmed by empirical data, is a very 
interesting issue for future research. 

A noticeable shortcoming of our analysis pertains the prediction of a linear impact with order 
size (Q = Th in our notations), which contradicts the observed concave behavior A ~ , with 
7 G [0.5,0.8]. It is easy to identify the origin of this behavior within the MG, in the assumed 
market mechanism. This assumes a matching of orders at a single deterministic price, fixed by 
market clearing. Indeed, the price itself does not enter into the definition of the MG. Strictly 
speaking the MG describes the dynamics of the excess demand, and how this is determined by 
market makers who adaptively respond to it, in their effort to match it dynamically. Therefore, 
on one side this suggest to look for alternative definitions of the price in the MG, capturing 
more faithfully the dynamics of the order book. On the other, this suggest to look at empirical 
measures of the market impact in terms of the demand imbalance itself. 

The analysis of the more realistic case in which the size of the orders hT is Pareto distributed 
would also be a relevant extension of the present work, and could allow the model to describe 
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persistence of the order flow (see Farmer et al. ( 20 111 ) and Toth. Lemperiere et al. ( 201ll )). There 
are two possible choices: either h can be Pareto distributed and T fixed, or, vice versa, T is Pareto 
distributed and h fixed. In the former case it would not be possible to apply linear response 
theory, since there will be realizations with /i ~ P, which is not a small perturbation. In the 
latter case it must be considered that when one measures the market impact averaged over orders 
of different volumes for a given value of volume Q, only the orders whose volume is larger than 
or equal to Q are taken into account. If the difference in volume depends only on the difference 
in the length of the order then, at a given time, only the orders which are still standing are 
taken into the average. As a consequence, it is easy to show that, at least during the execution 
of the meta-order, the market impact is still linear with respect to the volume. Alternatively, it 
would be very interesting to study a Pareto distributed order flow for the environment (i.e. by 
assigning to the producers strategies which lead to large fluctuations of the excess demand Ap). 
In this case the market stationary properties should be different, and an anomalous response to 
a meta-order could emerge. It is not easy to anticipate what the behavior of the system would 
be as the problem becon ies harder to solve analytically . 

Finally, as discussed in Toth. Lemperiere et al. ( 201ll ). the market impact is expected to cross- 
over to a linear behavior for times longer than few days. This supports the view that the MG 
provides a description of market activity on intermediate time-scales, with P representing a 
time-scale between one day and a week. 
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Appendix A: The stationary state of the grand-canonical Minority Game 
with replica symmetric approach 

The problem of calculating the macroscopic observables of the grand-canonical MG in the sta- 
tionary state can be solved by finding the minimum (also referred as ground state) of a Hamil- 
tonian function, whic h is precisely the Lyapunov functio n minimized by the dynamics of the 
grand-canonical MG ( De Martino and Marsili et a/.l 120061 ) 



= —Y,{A\^i)\ (Al) 

^ M=l 

where 

(^l/i) = ^/Wpx^ + <^imi (A2) 

i=l 

and m = (mi, 772-2, • • • , ^N,)- Each of the rrii € [0, 1] is the average value of (pi in the stationary 
state, and is a Gaussian variable with zero mean and unit variance. For convenience we are 
considering an Hamiltonian that is different from the predictability H = IhsHq. Therefore, in 
order to recover the free energy associated with H, one has simply to multiply the free energy 
obtained with Hq by 2ns. 
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Let US define the partition function 



Z = Tim exp(-/3i?o) = Tim ( exp j ^ i^^{A\fi) 



(A3) 



where we have used the Hubbard-Stratonovich trick in the second equality and 

1 



exp{-zy2)f{z)dz . 



(A4) 



The minimum of Hq averaged over the fixed strategies af can be obtained from the partition 
function in the limit of /? — )• oo. It is given by 



hm mm — 

N^^oo m A*, 



lim lim lim 



1 



■ In Z*^ , 



(A5) 



N.^oo /3-s-oo n-s-0 PNgU 

where the over-bar denotes an average over the fixed strategies and the variable n was introduced 
i n order to calculate the average InZ through the one of Z", a strategy known as replica trick 

(|Mezard et a/.lll987l ). We no w proceed calculating Z". 

It is a simple exercise (see De Martino and Marsili et al. 20061 ) to show that after averaging 
out the fixed strategies we obtain 



Z"" = Trjmj ( exp 



a,b \ i 



•miarriib 



(A6) 



where the indices a, h run over the n replicas. 

Introducing the overlap rnatrix Gah = 7^- Ylii '^ia^^jb and the response matrix Gah (see 
De Martino and Marsili et al. 20061 . for details), it follows that 



Z'^ = y dGciGexp(-/3nA,/(G,G)), 



(A7) 



where 



/(G,G) 



-Tr In 



InTrr I exp ( ^ Gabmafrib 



(A8) 



We now take the ansatz (known as the replica symmetric ansatz) Gab = 9 + {G — g)5ab and 
Gab = 2i? — {R + p/(3)6ab- Note that in terms of the original variables 



1 ^= 



(A9) 



i=l 
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In this way the free energy /(G, G) becomes 

-^(inj^ dmexp{-(3V{m,z))^ , (AlO) 

where y(m, z) = ^;-pn^^ + \J^^f^- the Umit of — )■ oo the free energy is dominated by its 
minimum and the saddle-point equations are 

G = F(C) 



1 + X 2 
1 

P 



^Erf(C/V2) 



1 + X 

where x = PiG-g), C = p/^/^R^, and F(C) = exp(-CV2)+ 2J^Erf(C/^/2) + iErfc(C/^/2). 

Rearranging these equations we obtain the following recipe to calculate the macroscopic observ- 
ables. First one should solve the following transcendental equation in order to obtain as a 
function of the control parameters and np, 



n,F(C) = ^ - . (A12) 



After that, one can express the macroscopic observables as function of C, by using the saddle 
point equations. For example, it is easy to calculate the minimum of the free energy to obtain 
the predictability as 



H = 2nlfsP P = Ns {np + 7i,F(C)) ( 1 - ^Erf(C/\/2) ) ' , (A13) 



where fsp is the free energy at the saddle point. Moreover the susceptibility x is written as 

X = ^^iM£A^. (A14) 
2 - „,Ert(C/V2) 



Appendix B: Perturbing the Hamiltonian 

We now consider the perturbed Hamiltonian 

p 

2N' 



HH{m) = —Y,{{A\p) + h)\ (Bl) 



/.=i 



22, 2012 4;14 Quantitative Finance paper 



14 



This corresponds to the Lyapunov function that is minimized by the dynamics when the meta- 
order is turned-on. Therefore, it gives the stationary state measure for the case with an extra 
buyer. 

The calculations in the present case are analogous to the calculations in Appendix El For 
example, with the perturbation, Eq. (IA6j) changes to 



exp I 



^ " a,ij. / a,b \ * i 



, rriiamib 



, (B2) 



where the Zh indicates the partition function associated with Hh- This leads to the following 
free energy 



(B3) 



a, 6 



where Mab = Sab + P{np/ns + Gab) and /(G,G) is the free energy given by (jAlOp . With the 
replica symmetric ansatz we have ^ab-^ab ~ i+p(G-g) + leading to 



h^P 



+ f{G,g,R,p) 



2iV|(l + x) 

Hence, we have obtained the free energy for the perturbed case. Moreover, if we write 



(B4) 



= Tr™ exp ( J] {{A\fi) + h) 



(B5) 



we get 



lim lim 



This gives 



d 



(B7) 



Finally, from Eq. ()B4p we obtain Eq. (|12|) . where x is given by Eq. (|A14p . after numerically 
solving Eq. (jll2]) . 



Appendix C: Permanent impact 



In the following we show how the permanent impact relates to the stationary state quantities of 
the model calculated above. We will consider the continuous time description of the dynamics 
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(IChallet et fliJl2nn,^ ) as justified in Sec. [Sj and define accordingly the quantities T = T/P and 
A(t) = A(tP) « j^^ dr' {a,,{t'P)) - (A;,=o(t'P)) • (CI) 



Within linear response theory, Eq. (|13p yield the response at time t of the system due to a 
perturbation /i(t), in terms of the kernel K{t). Note that J^^cLtK^t) = 1 and that causality 

requires K{t) = for r < 0. From Eq. Q the permanent impact A* = limT-_j>oo A(r) is given 
by 



oo 



A* = I dT I dr' K{t - t') h{T') . (C2) 



oo 



Assuming that K{t) is an integrable function, we have 



-+00 J 



Therefore, the expression for the permanent impact can be written as 



h , , , /■+°° du 



dr / dr' / _K,,e-^'^(^-^') 



(1 + X) Jo Jo J~oo 27r 

h f+°° duj / 1 \ /l-e(*'^+^)^\ , , 



(1 + x) <s-s>o+ y_oo 27r \iuj + e J \ —ioj — e 



The value of this integral is determined by the value of the residues of the integrand on the 
plane Quj > 0, provided that \Ki^/uj\ — t- for \oj\ — t- oo. If K^^ has no poles on such domain, 
which corresponds to the assumption of causality K{t) = for r < 0, the only contribution to 
the integral comes from a second-order pole in uj = ie. Hence, using J_ °° dTK[T) = 1 we obtain 
our final formula 

A* = i-^ hm (-iK,,f) = . . (C5) 



Appendix D: Average execution cost 

We calculate in the following the expression for the average execution cost as defined in Eq. (Ilip , 
and compare it with the maximum price increment A(T)p(0). To calculate those quantities it 
is sufficient to assume linear response by plugging Eq. (|13p into the definitions of A(T) and A. 
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The former quantity can be obtained straightforwardly integrating by parts, and results 

MT)^J^ r dr r dT'Kir-T') 
1 + X Jo Jq 



h 



1 + X 
1 + X 



T 



dT (T - r) K{t) 



(Dl) 



where the averages have been defined in Eq. (jl6p . The expression of the average execution 
cost is analogously derived by using integration by parts 



A = l A{t)^l TdrA 

l<t<T 



(r) 



2(1 + X) 



12 



(D2) 



By taking the ratio of those expressions one obtains Eq. p6p . The functions can formally 
be computed by inserting the expression (jl4|) into their definition. In particular we obtain the 
series expansion 



— ' n\ n + m + 1 

71=0 



(D3) 



where Ki^\^) is the n-th derivative of the function Kj.{t) evaluated in zero. Such expansion 
allows to calculate Eq. (I17p which is obtained in the limit of small The execution cost in the 
limit of T much larger than the characteristic relaxation time can be estimated by noting the 
scaling 



1 /r, 



m+1 



(D4) 



so that just the terms contributes to Eq. (jl6p . This implies that in this limit A/A(T) tends 
to the value of 1/2. 
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